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A Method of Developing the Perturbative Function of 

Planetary Motion. 



By Simon Newcomb. 



The development of the perturbative function in powers of the eccentricities 
and inclinations has of late been generally regarded as of little value, owing 
to the complex character of the series to which it leads. It is, in consequence, 
but little used, even in those cases of nearly circular orbits where it would be 
most convenient. Still, it is the only development in which the disturbing 
force is given as an explicit function of all the elements, and is therefore of 
more interest to the geometer than any other. Moreover, it admits of various 
simplifications in its application to the numerical problems of celestial mechanics 
which deserve more attention than they have received, and which may entitle 
it to a more favorable comparison with other methods than it has been supposed 
to offer. 

The object of the present paper is to exhibit a method of effecting the 
development in powers of the eccentricities, which seems to me to offer some 
features of interest, and possibly to contain the germ of some principle which I 
have not fully grasped, and which may admit of wider and more important 
applications. I refer especially to the expression of the coefficient of each 
power of the eccentricity in terms of the coefficients of lower powers, and to 
the expression of the coefficient of each term involving the perihelia of two 
planets as the symbolic product of coefficients involving the perihelion of one 
only. The first of these features was pointed out in a note to the French 
Academy, found in the Comptes Rendus, Vol. LXX. p. 385, the ground of which 
is covered by the present paper. The second was discovered on the completion 
of the theory some years later. 

One great practical advantage of the process is that it is reduced to a 
uniform operation of algebraic multiplication, which can be executed by an 
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unskilled computer, and can be carried to any extent without repeating the 
previous processes. I see no reason why it might not be possible to express 
the terms of the moon's longitude by some similar series of operations, and 
thus greatly simplify the practical problem of the lunar theory. Probably each 
step would be found to involve the solution of a differential equation, but this 
equation might be of a very simple character, and only a particular integral 
would be required. 

The present development offers nothing new in the method of taking 
account of the mutual inclination of the orbits. The development with respect 
to this element may be made by any method which gives the coefficients as 
known functions of the radii vectores and the mutual inclination, and the ele- 
ments under the signs sine and cosine are the multiples of the distances of 
the two planets from the common node. 

The expression to be developed is 

T 

R = (r 2 — 2rr' cos V — r' 2 )~ 4 — ^75 cos V, 

in which r and r' are the radii vectores of the planets, and V the angle between 
these radii vectores. 

The second term of this expression admits, after the development is effected, 
of being merged in the first by a simple and well-known modification of certain 
terms of the first ; we shall therefore confine our attention to the first term. The 
following notation is used : 

V, v', the true angular distances of the planets from their common node. 
A, A.', the mean values of V and v'. 

y, the mutual inclination of the orbits. 

o-, sine \y. 
p, p, the logarithms of the radii vectores. 
», »', the logarithms of the mean distances. 
e, e', the eccentricities. 
g, g', the mean anomalies. 

a, the ratio of the mean distances. 

By substituting for cos V its known value, 

cos V= cos V cos V' + sin V sin V cos y 

or cos V= (1 — ex 2 ) cos (v — V') + o- 2 cos (v + v'), 

and then developing in cosines of multiples of V and V', we shall obtain the 
terms with which we are to set out. In order not to weary the reader with 
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what is not essential to the present object, we shall state only those conclusions 
which form the basis of the new method. 

If we suppose the eccentricities to vanish, the value of R can be developed 

in the form 

R = \ 2,Ai cos (iX' — iX) 

+ ZB t cos ((» + 1) X'- (i - 1) X) 

+ 2a cos ((» + 2) X' - (i - 2) X) (a) 

+ etc. 

where the index i takes all integral values from positive to negative infinity, 
and Ai , B { , C { , etc. are functions of the mean distances and inclinations which 
admit of explicit development in powers of cr. 

This expression for R may be thrown into the form 

2? = 2,2,A„ 1 COs(vX'+ /t X) (1) 

where ju, and v each assume all values from — co to + co , but not independently, 
being subject to the single restriction that both values must be even, or both odd. 
The coefficients A are homogeneous and of the degree — 1 in a and a', 
admitting of being expressed indifferently in either of the forms 



5* (I) ° r a*(f)' 



It is in practice more convenient to choose the form in which the fraction —, or ~ 

shall be less than unity, but, for the purposes of the present investigation, the 
choice is indifferent. 

The original expression for R is of the form 

R=f{v,v',r,r',a); (2) 

and, assuming the eccentricities of both planets to vanish, we have supposed it 

developed in the form 

R =f(X,X , ,a,af,a-). (3) 

In order to express R as a function of the eccentricities and other elements of 
the two planets, we must now substitute V, v', r, and r' for X, X', a and a' in (3), 
or, which is the same thing, in (1), the quantities to be substituted being ex- 
pressed in terms of the eccentricities and mean anomalies. 

To continue the process we shall consider R as a function of the logarithms 
of the radii vectores of the planets, instead of the radii vectores themselves, 

putting 

p = log r, p = log /, 

r> = log a, »'= log a', (4) 
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. da dr 

which gives ^ = a ; Tp = r. 

It is not necessary to have either A or its derivatives expressed explicitly as a 
function of *> and »', since we have, with respect to any function <£ of a and a', 

dcj> d<f> . dj> , d<j> 

dv~ a da' fa'~ a da'' 

a form hy which -all the derivatives with respect to » and »' may he expressed as 

functions of a and a'. 

All the coefficients of A^^ heing homogeneous and of the degree — 1 in a 

and a', we have 

— _i_ — 4- a — n 
d-o "r dv' "T" -^ ~~ u - 

We shall hereafter use a symbolic notation, putting 2? for the operation ^ and 
JO for ^7 . All the derivatives of A with respect to » and »' being, like J. itself, 
homogeneous and of the degree — 1, we may put, in general, 

D + D' = -l, (5) 

and may combine D and D' as if they were multipliers according to the usual 
rules for such symbols. 

To effect the development we require we must, in (2), put 

V=\ + <$>{e,g), 
p =r> + $(e,g), 

F'=X'+<He',<7'), 

p'=»'+W,sO; (6) 

the function <f> representing the equation of the centre, and i// the portion of 
log r which depends on the eccentricity. From these equations we see that 
considering R first as a function of V,v',p, and p, and then as a function of 
X,'X',» ,r>',e,e',g and g',we shall have for any compound derivative with respect 
to v and p 

Jm + ml + n + it' r> ^m + m' + n + n'JJ 



d V m d V' m 'dp n dp m ' d\ m d\" n 'd-D n d-D' a ' 

that is, any derivative with respect to V ,v', p,or p is found by taking the corresponding 
derivative of the developed function with respect to X , A.', » , and »' 

Now, supposing that, in (2), V and p are replaced by their values in (6), we 
shall have 



Newcomb : On the Perturbative Function of Planetary Motion. 197 

dB dB dv . dB dp 

de dv de dp de 

dv dB . dp dB /k\ 

de dk * de d?> ' 

This equation is the fundamental one in our method of development. By it we 
express the derivative of R with respect to e in terms of its derivatives with 
respect to A. and v. Let us now differentiate this expression with respect to e 

n times in succession, representing by _D A the operation ^. We thus obtain 

d n+1 B _ ( dv d^B /n\ d?V d n - x B /n\ d?V d n ~ 2 B ) 

de n + 1 ~^ \de de" "*" \1J de 2 de"- 1 "+" \2J de* de n ~ 2 "•" \ 



L \dp d n B 

+ ^4, \fo -^ 



fn\ d*p dr-JB fn\ d?p d n ~ 2 B ) „ 

+ \1J de 2 de*- 1 + \2J de s de"-* + \ W 



Thus, we have expressed the derivative of any order with respect to e in terms 
of the derivatives of lower orders, and, by successive substitutions, this derivative 
will be expressed in terms of derivatives of R with respect to X and r>. 

The coefficient of e n+1 in R is found by putting e= in (8), and dividing 
by 1 . 2 . 3 . . . . n + 1. In strictness we should suppose e = only after differen- 
tiating with respect to A, but it is evident that the two operations may be inter- 
changed without affecting the result. If we represent this coefficient by R n+1) , 
and replace the derivatives with respect to e by the corresponding values of this 
coefficient, namely, 



dB 
de 



Rv 



^=21B» 



l£ = 3!22« 



etc. etc. 

If, also, we represent, as in these last equations, by the subscript zero the 
operation of putting e = after differentiation, we shall have 
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2g<»+i>= n 



j n\ dvo R(n) , /n\ (»-!)! dV„ 

| (« + l) ! de M ^ ^lj + 1) ! de 2 •** ^ etC ' 

^■^-Wn+l)! de -^ ^ll J (» + l)! de 2 -« 1" etc. 



If we represent by v n and p n the coefficients of e n in F and p respectively, we 
shall have 

^ e n n. y n ; de n re . /)„. 

Substituting these values in the above equations, we find by simple reductions 

(n + 1) B-+ 1 = D K { vJBP + 2 Fafle -« + 3 v s W n ~ 2 > + .... + (» + 1) r n+1 E } 

+ Z>„ { Pl ^ + 2p 2 B n -v + 3p 3 R( n -v + .... + (n+l)p n+1 B }. (9) 

Thus, each coefficient is expressed as a linear function of the derivatives of the 
coefficients of lower orders. We have next to substitute for F 4 and p t their 
values in terms of g. We have in general 

*F, = *T Ttpsinjg; i Pi = hf cob jg (10) 

from the developments of the elliptic motion, but it is to be remarked that j 
does not assume all values between the limits + i and — i, but only every alter- 
nate value. The special values of k and h to terms of the seventh order are 
shown in the following scheme. The values for negative values of j are formed 
from those for positive values by the formulae 

h% , = - kf ; Wj = hf. 

In strictness it is not necessary to suppose j negative at all, since the complete 
value of 4 p 4 , for example, is 

4/> 4 = iK Y + hf cos 2^ + hf cos ig, 

but it will be found convenient in forming the required functions to suppose j 
negative in some cases. 
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Values of h and k. 



+ 2 


-l 


k'-i 
-2 


-1 














5 

+ 2 


h\ 

3 

2 





K 
+ 1 


^—2 

5 

~~2 


A'l 2 

3 

~ 2 








13 

+ T 


A3' 

17- 
8 


3 

4 


AT 
9 

+ 8 


3 
+ 4 


9 

+ 8 


k'l s 
13 

~ 4 


A_ 3 

17 
_ 8 


k\ v 

103 

+ 24~ 


AJ V 

71 

24 


11 
"~~6~ 


hV 
11 


7,1V 




AJ V 

1 
+ 4 


Z.IV 
11 


Ai v 2 
11 

+ T 


1097 
+ 192 


523 
~128 


kj 

215 

~~ 64 


A 3 V 

385 
+ 128 


kj 

25 
+ 96 


Af 

5 
+ 64 


kit 

25 
96 


5 
+ 64 


kj 1 

1223 
+ 160 


hj 1 

899 
~160 


kj 1 
451 


hj 1 

387 
+ 80 


17 

+ 32 


9 

~32 


A„ VI 



A VI 

1 

+ 8 


A,™ 

330911 
+ 32256 


hj 11 

355081 

46080 


kj* 

41699 
4608 


hj* 

70273 
+ 9216 

1 


* 3 vn 
665 
+ 512 


A3 711 
5201 
5120 


kj n 

749 
+ 4608 


Ar 

889 
+ 9216 



To commence the development, let us take any one term of (1), 
brevity put 

and omit writing the indices /x and v. The term will then be 

B = B m = AcosN, 
and its derivatives with respect to X and » will be 



and for 



cfa> 



= 2X4. cos iV. 
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The general formula (9) gives, by putting n = 0, 

^ = F i ~dx + Pi sr 

= — fik[A sin g sin iV+ h[DA cosg cos JV, 
whence 

2 E a) = (p%A + h[DA) cos (iV + y) 

+ (— ^i + AJI>il) cos(i^— gr). 

By the repeated application of (9), supposing w successively equal to 1, 2, 3, etc., 

we shall obtain the successive values of R n) . This process will be facilitated by 

finding a general formula for passing from R {n) to P in + 1) . For this purpose let 

us put 

Pj, the coefficient of cos (iV + jg) in E {n) , 

we then have, in the special cases n ' = and n = 1, 

PJ = (/**£ + h[D)A, 

FLi = (-fi% + KD)A, 
and in the general case 
jgw = p» cog ^ + nflr ^ _|_ p»_ 2 cog (_^r + (n _ 2 ) g ) + . . . . + Pl„ cos (iV^— n#), 

the index j taking each alternate value from + n to — n. Differentiating with 
respect to X and », we shall have 

D K R n) = — /*P£ sin (iV+ tip) — /*PS_2 sin (iV + (n — 2) g) — etc. 

2>.i? w = DPI cos(i\T+ ngr) + X>P^_ 2 cos (i\T+ 0» — 2) g) + etc. 

Putting side by side, and writing in the most condensed form, the pairs of factors 
which enter into (9), we find them as follows : 

, j = + n 

D x R^ = -fx, 2 PJ sm {N+jg) V 1 = k[smg 

j = —n 
j = n — l 

D k R a - l > = - l i 2 PJ- 1 s\n(N + jg) 2 v 2 = 1% sin 2 g 

j=—n+l 
j = « — 2 

D x R n -^ = -fj, 2 P;- 2 sin(iY r +,;V) 3 F 3 = #" sin 3 ^ + #" sin # 

i — — n -f 2 



i = n -f-1 

Z> A 2?°> = — /ti'gsiniV" (» + l)r, +1 = |2 $' 1 + a) sin ^ 

i = — » — 1 
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j = n 

jDhpjg.*) = 2? 2 PJ cos {N + jg) p 1 = h' 1 cosg 

j = — n 
,?' = »— 1 

jtoft—v = D 2 PJ- 1 cos {N + jg) 2 p 2 = hi cos 2 g + J %' 

DoBS*-* = D~t . PJ' 2 cos {N + jg) Bp 8 = 1Q cos 3 # + J£ cos ^ 

j = — «-f-2 



X>»i2(°) = Z?J*> cos N (n + 1) p ra + x = * |T ft? + a) cos ty. 

t= — » — i 

If we form the products of these quantities according to the formula (9), and 
compare the coefficients of the several angles, we find 

2(n+l)P:tl = (pK+h' 1 D)PZ 

+ (pic; + kb) pi-\ 

+ {^+JhD)Pl-\. 
+ 

+ {^mtv + Ki%yD) p% 

2 (n + 1) Pit \ = (/**; + ft-D) PS-2 + (M-i + ALiD) -PS 

+ (^ + Ai'i?) pszi + ^'i>p:=1 

+ etc. + W + £"!>) P:=5 (11) 

+ (/*%•> + AS ,) 2?)i^.i +etc 

+ (hW±? + v:±\w)pi 

2(n + l)i*±J = (/J54 +^)i*U+(/**Li+AliD)i 3 S_, 

+ (/*# + cd) pr^+ kd^z\+ {tiicu + ^- 2 -0)'pri 

+ W + cc>) Pri+ W + Ai"i>) Pti + (^-i + ki x d)p::i 



To show the law of progression of the several classes of terms, I have purposely- 
written some of the coefficients k and h with negative indices instead of using 
the corresponding positive ones, which, in the actual development, will be 
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substituted for them. The law can be seen by induction and comparison with- 
out a full statement of it. As we diminish the value of the lower exponent of 
P by successive steps, two units at a time, the number of columns of products 
increases by one at each step, but the number of products in each column 
diminishes in consequence of one of the factors vanishing whenever a lower 
exponent in h, Jc, or P exceeds the upper one in absolute value. Supposing the 
successive values of PJ + 1 to be written in this form with continually diminishing 
values of j until we reach the value j = — (n -+- 1), to this value would corre- 
spond n + 2 columns, of which, however, the first would entirely vanish, as 
would every term of the remaining ones, except the first. Placing these terms 
in a column, they will be as follows : — 

+ (/**", + h'L^Pl-^ + (- txTcl + KB)P^-\^ 

+ (M'-s + h'l 3 D) Pl- 2 +2 + (- f,ks + %'D) Pl- 2 +2 

+ + 



This is a particular case of the general law of formation of the values of P when 
j is negative, which law may be -expressed as follows: — 

Each value of P with j negative may he formed from the corresponding value for j 
positive by changing the sign of p, and of the lower exponents in all the values of P which 
enter into it. 

We have next to consider the development with respect to the powers of e'. 
What we have hitherto done has been to take 

R=f{v,p,v', P ') 
and substituting 

F=X+^(e,gr)j 

p = » + i|», (e , g), 

to find the coefficients R w of the development with respect to e only. We have 
written X' and »' instead of V and p in giving the values of the coefficients. 
Now, when we replace X' and »' by v' and p , expressed in terms of the elements, 
each coefficient JS (n) will become a function of v' and p, and hence of e', g', etc. 
We shall now proceed to develop each of these coefficients in powers of e' by the 
same process which was followed in developing R in powers of e. Let 
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be any term of B {n) , and, having replaced X' by v' and »' by p in this term, let 

us suppose 

F ' = X'+<£(e',/), 

p' =* + +(<{,&; 
and let us put 

i? 1 '"', the coefficient of e v in this development 

We shall then have, as in (9), 

(»' + 1) ^•»'+ 1 = v{D y E n ' n ' + 2 r^DyPJ 1 ^- 1 + ....+ (n' + 1) V^ +1 R n '° 

+ p' 1 JD'B n ' n/ + 2 p^'i^'- 1 + .... + (»'+ 1) P ' n/+1 R tt '°. (12) 

The expressions for f' and p' will be the same as those for F and p in (10), except 
that </ is to be substituted for g. We may now investigate the general law of 
development in the same way as before. Putting, for brevity, 

N' = pX + v\' + jg = N+jg, 
so that 

P; cos N' 

is any term of R in) or, which is the same thing, of M n '°, let us represent the cor- 
responding terms in P n ' n ' by 

£n,n, =p n,% cos ( N ' +n Y)+P« : *_s COS {JT+ K~2)/) +. . . .+P];%ZOs{N'-rig'). 

Then, proceeding as before, we shall find 

2 (n' + 1) P];^U = (v% + KIT) P" t % 

+ (v% + %'&) P*j#z\ 

+ 

+ (vWi\> + w$iw)P n j: i (13) 

All the other coefficients can be formed from the corresponding ones of (11) by 

writing 

v for p, 

D' for JD, 

P n j for P, 

n' for n. 

We have now an important law of this development of R to bring out. First, in 
the equations (11), by supposing in succession n= 1, n = 2, n= 3, etc., and by 
continually substituting in each set of terms the values of those of a lower order, 
we shall finally express all the values of P 1 ) in terms of P\ = A and its successive 



204 Newcomb : On the JPerturbative Function of Planetary Motion. 

derivatives with respect to •». Moreover, the operation of forming these deriva- 
tives being always linear, we can combine all the operations represented by the 
symbols /xk + hD as if D represented a coefficient ; that is, having the quantity 

(pk + hD) (pi? + h'D) ....(fik n + h n D) i*> 

we can multiply these several symbols as if D were a coefficient. By this opera- 
tion, putting A for P%, we shall finally obtain an expression of the form 

P] = WjA, 

in which II" represents an entire function of p and D of the degree n. 

Secondly, by treating the equations represented by (13) in the same manner, 
we shall be able to represent each value of P" ; y in the form 

rm n' — TT"' ( T> n > — jm\ 

r 5 , y 1 V K-^j > o * 5) 
in which 11^ represents an entire function of v and D of the degree ri. Substi- 
tuting for P'j its value just given, we shall have 

P>j$ = H$IFjA (14) 

in which the two symbols can be combined by the rule of multiplication. It thus 
appears that when we have found the development in powers of e for q' = 0, and that in e' 
for e = 0, we have solved the whole problem, and the terms multiplied by any product of a 
power of e by a power of e' can then be found by a symbolic multiplication. 

It may not be amiss to recapitulate the result which we have reached. Sup- 
pose that we develop R in powers of e on the supposition e' = 0, and that any 
term of this development is represented in the form 

R = e n ILjA cos (JST+jg), 

A being a function of the mean distances, II", an operating symbol, and JSF a 

function of X and X' which does not contain g. Suppose, next, that we develop 

R in powers of e', putting e = 0, and that any term of this development is 

represented by 

R = e' n 'n f 'A cos (JST + j'g'), 

then the coefficient of cos {N+jg + j'g') in the complete development will be 

represented by 

e n e' n 'UJU f 'A. (15) 

To proceed to the actual development, it is necessary to form the symbolic 
factors represented by 11" and 11"/. This we may do from the forms (11) by 
substituting for h and k their numerical values, and substituting the symbol 
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nj for Pj, remembering that n® = 1. We thus have, for each successive 
value of n, 

2(n + i)n:tI = (2/*-2>)ns 



+ (§/*-§tf)n:zl 

, /13 17 ^\ 

, /103 71 ^\ 

/. 1097 __ 523 \ tc _ 4 
+ (^192 /* 128 D j^ n n -4 

/1223 __899 \ „_ 5 
+ (^160 ^ 160-^J n «-5 



+ etc. 
the series terminating with Hq. 



2(n+l)n£i = (2 / t-z>)n5_ s +(-2 /t -i>)n: 
+ (§/*- 1 2>)n:=i + i>n»zj 

+ (^ - £u) n-zj + (- ^ + §z>) n-zS 

+ etc. + etc. 

+ (^V + A«2>) nz a + (#l±¥ii + h*±¥D) III 

2{n + l)Til±\ = (2p-D)Tll_ i +{-2 l L-D)W n _ 2 

+ (^-§z>)n^ + i>n»zi + (-|^-|z>)n»zl 

+ (^-yD)nn +(-!^+^)n»zf + (|^ + |i))n:zi 

+ etc. + etc. + etc. 

+ (»/* + h^P)uU + (*SV + A$V>)ni! + (i^s+S'/t + A ( ^Z))ng 
2(» + i)n^ = (2/* - i>)n»_ 6 + (-2/* - 2>K_ 4 

+ (|/t-|/>)mzj + i?n»zJ+(-|^-|i>)n»z^ 



2 

w — 4 



+ 



(-^-^d)ii»: 



+ etc. + etc. + etc. 

+ (J£zf/*+ A<fc#2>)n!U + (^zi>+A'^'i>)n 2 _ 2 + (*SV+A£V?)iHi 

+ (*:+?/* + A«:±s>i?)ii8 
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2 (n + 1) iri+Li = (- 2'/*" - D) rc. 

+ (-l/--l i> ) n -^ 1 +l. 

+ (-T^-T^) n -+ 2 

+ etc. 

+ (-#:#/*+ Wi)) OS 

2(w + l)n!LtVi = (-2 i a-i>)ni K + 2 +(2 / ,-D)ni„ 

+ (-^-|x>)ni- 1 + a + .diel» 1 + i 



— 2 

+ 2 



+ etc. + etc. 

+ (-»?/* + #»>D) nj + (- ui±¥p + tt:±yD) n° 

The actual computation of II" for negative values of j is not necessary, since its 
values may be obtained from those for positive j by changing the sign of /jl. 
These numerical coefficients may be continued to any extent by means of the 
scheme of values of h and k already given, it being remarked that the successive 
columns appear in the same order as in the scheme. 

We are now ready to proceed with the actual computation of the functions 
P n j ; % , or, which is the same thing, of the symbolic functions IT" ; % which express 
the values of P} ; y when considered as operators on the functions A { , B t , etc. 
We shall give only a few of these functions for the purpose of illustrating the 
method. 

We consider, firstly, the general values of II" ; <J and of IIq ; % which arise when, 
in the general term of R for circular orbits, 

R = A cos (v\' + /uA), 

we substitute for the mean values of the radii vectores and longitudes those which 
correspond to the elliptic motion. The functions of II for values of j and / to 
the fourth order, inclusive, are then found to be as follows. It is not necessary 
to write the functions for negative values of j or /, because they are found from 
the corresponding ones for positive values by simply changing the sign of /a and 
D, but a few are given for perspicuity. 
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rig =1 



UU = -p.-oD 



n 2 =-^ 2 +| d + \& 

n 2 _ 2 = ^-gV + (gft-^D + gD" 

HLt = g/t 8 ; - §V + g> + (^ 2 - jg/* + jg) -D + (- g> + ^) D 2 - ^D 3 . 

TT4 ! 4 I 5 3 _L 283 2 _1_ 103 I / 1 3 ! 2 5l 71 N 71 

n I =24/*+i6/ t +384/ 4 + 192^+ ^~~ 12 /* 2 /* — 64 /* — 192J -° 

ns =i^-M^ 2 + i^+(-^ 2 -6i) I)2 -i^ + 6i^ 

The values of II' being formed, as already shown, by simply changing //, into v 
and D into D', it is unnecessary to write them as functions of D'. As it will be 
convenient to have but one form of derivative, they should be transformed, the 
symbol 2/ being replaced by D by writing 

jy=-(l + D). 
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We thus find 



II°;i = v + | + |z> 



nS:5 = ^ + |i' + 2 + (sv+-g)D + gD» 



nH = -^ +\d + \& 



U *:* — U v + 48 ^ + 384 ^ + 192 J ' + 3 + V 12 8 64 6 V 

llo.a— e^ 24 v 6^ 48^ & \ & 2 V 96 v 48 J n 

The symbolic products 

are functions of both the indices /* and p. Instead of using them in their general 
form, it is more convenient to apply them to the separate terms of the develop- 
ment (a). To develop the first term, we put 

fji = — i ; v = i. 

For the second, 

fj, = — i+ljj/ = *+l, etc. 

To show the forms to which we are thus led, the following exhibit of the first 
four orders of terms arising from the first term of (a) is presented. The symbol 
(i) indicates that in II, fi is changed to — i and v to i. 
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Developed Value of P. 
B= i {1 + e 2 ng + e'X; 2 , + ^n 4 + eV*n8;8 + e' 4 ng; 4 + etc.}' ' A g 

t = oo 

+ S 

i = l 

{1 + e 2 n 2 + e' 2 ng;? + e 4 nj;J + eV 2 n 2 ;§ + e' 4 ng; 4 + etc.}"' J, cos («V - iX) 



+ e {n} 
+ e' {ng 

+ e 2 {nr 

+ e' 2 {ng 

+ e 3 (n? 
+ eV {n 2 
+ ee' 2 {n} 

+ e 3 {ng 

+ e 4 {n 4 
+ eV {n 3 



g+e'tt; 



+ eV 2 {n| 
+ ee' 3 {n} 

+ e' 4 {ng 



S + e 2 n 3 
i + e 2 ng 

g + e 2 n| 
} + e 2 n? 

l + e 2 ni 



1 = -f oo 

+ a 

* = 00 

g + e' 2 n{ 
J + e' 2 ng 



i + e 2 n 4 

l + e 2 ng 



2 + e 2 ns 



1+ e 2 nf 
l+e 2 m 



g + e' 2 ni 
i + e' 2 ni 

I+e' 2 ng 

S+e' 2 n§ 
I + «*n3 

I + e' 2 nj 

3 + e ,2 ng 
g + e /2 n| 

J + e'^l 

l + e' 2 ni 
l + e^nj 
t + ^ng 



§+e 4 n?;g+eV 2 n?;i+e' 4 ni 
!H-e 4 n 4 ;J+eV 2 ng;f + e' 4 ng 

I + *nj : 8 +^v«n};8 + e' 4 ni 

f+e 4 nf;i+eV 2 n?;?+e' 4 nj 
|+e 4 n 4 ;I + eV 2 ng;| + e'X 

I+e 4 n?;g + eV 2 nf;„ 2 + e' 4 n! 
3 +e 4 n«;l + eV 2 n|; 8 +e' 4 n 2 
| + e 4 n?;l + eV 2 nf;|+ e' 4 ni 
|+e 4 n 4 ;l + eV 2 ng;|+e' 4 ng 



4 + etc.} (i) Af cos O'A' — iX + g) 
f + etc.} (i) A ( cos (iX' — iX + /) 

$ + etc.} (i) ^ cos (iX' — i\ + 2 g) 
{ +■ etc.}'*' A t cos (*X' — t'X + g' + g) 
jj + etc.}'*' ^ cos (*V - iX + 2 #') 

J + etc.}"'' J, cos O'X' — iX + 3 g) 
I + etc.}""' At cos (*X' - i\ + g' + 2g) 
% + etc.} (i) 4' cos (i\' — i\+2g' + g) 
I + etc.}"' A { cos (iX' - *X + 3 g 1 ) 



I + etc.}"' At cos (iX' — ik + 4 gr) 
f + etc.}"' At cos (tX' — i\+g' + 3g) 
t + etc.}"' 4 cos (iX' -i\ + 2g' + 2g) 
5 S + etc.} "> ^ cos (*X' - *X + 3 g' + g) 
\ + etc.}"' 1 At cos (tX' — iX + 4 g-') 



This development is not directly comparable with those hitherto executed, 
because we use the successive logarithmic derivatives D, instead of the deriva- 
tives with respect to a, the ratio of the mean distances. The two classes of 
derivatives are, however, connected by a simple linear relation which makes it 
easy to pass from one to the other. 



